2D quantum gravity from quantum entanglement 
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In quantum systems with many degrees of freedom the rephca method is a useful tool to study the 
entanglement of arbitrary spatial regions. We apply it in a way which allows them to back-react. As 
a consequence, they become dynamical subsystems whose position, form and extension is determined 
by their interaction with the whole system. We analyze in particular quantum spin chains described 
at criticality by a conformal field theory (CFT). Its coupling to the Gibbs' ensemble of all possible 
subsystems is relevant and drives the system into a new fixed point which is argued to be that of the 
2D quantum gravity coupled to this system. Numerical experiments on the critical Ising model show 
that the new critical exponents agree with those predicted by the formula of Knizhnik, Polyakov 
and Zamolodchikov. 
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Understanding the effects of quantum entanglement of 
systems with many degrees of freedom such as quantum 
spin chains or quantum field theories is a challenging 
problem which connects statistical mechanics to quan- 
tum information science. If the system under study is 
taken to be in a pure state \^), a complete description 
of the information available to an observer who has ac- 
cess only to a subsystem A will be given by the reduced 
density matrix 



PA trs |*)(*| 



(1) 



obtained by tracing over the degrees of freedom of the 
remainder B, inaccessible to the observer. 

There are various functions of pA which could be used 
as useful probes to measure how closely entangled, or how 
'quantum', a given state is (see [1] for reviews). Most 
of them are expressed in terms of trp^. For instance, 
the Renyi entropy is RA{n) = logtrp^/(l — n) and the 
entanglement entropy is the limit Sa = liiii„_>.i RA(n) = 
-tr PA log PA- 

In quantum field theory the quantity tr for inte- 
ger n can be computed without explicit knowledge of 
the ground state through the so-called replica method 
0-13] ■ Following this procedure the partition function Z 
of a d dimensional quantum system is computed in the 
standard way by doing the functional integration on n 
copies (or replicas) of the corresponding Euclidean clas- 
sical system in d -I- 1 dimensions. These copies interact 
among themselves through the inaccessible subsystems 
B as explained below. In this set up the quantity tr 
is proportional to the canonical partition function Z„ (A) 
of the coupled system of n replicas. More precisely we 
have trp^ = Z„(A)/Z". 

In most previous studies the accessible subsystem A is 
chosen to be fixed and a major goal is to investigate how 
tr p\ and the entanglement entropy depend on A. The 
point of view which is taken in this paper is different. 
We treat A (or equivalently B) as a back-reacting, dy- 
namical, subsystem whose position, form and extension is 



determined by its interaction with the whole system. We 
implement it by 'summing over histories', i.e. by putting 
the system in equilibrium with the Gibbs' ensemble {A} 
of all possible subsystems. 

As we will see, when the system is put on a lattice the 
sum over the ensemble {A} is unambiguously defined. 
An interesting consequence is that the sum over {A} de- 
fines a non-trivial, translation invariant, modification of 
the system uniquely generated by the quantum entan- 
glement of the set of all the accessible subsystems. This 
is particularly interesting when the system undergoes a 
critical transition. 

One important question to address is whether the in- 
teraction with the ensemble {A} changes the universality 
class of the critical system i.e., in renormalization group 
language, whether the coupling to {A} is relevant when 
the critical system is in a pure state. In this paper we 
answer this question for the class of 1 + 1 dimensional 
quantum systems described at criticality by a relativis- 
tic, massless, field theory, i.e. a conformal field theory 
(CFT) with central charge c. It turns out that any fixed 
point with c > becomes unstable when coupled to {A} 
and that the system flows to a fixed point with c — 0. A 
series of numerical experiments on the Ising model indi- 
cates that for a suitable choice of the coupling parameters 
the new fixed point is critical. The scaling dimensions of 
the primary fields turn out to be those expected in the 
coupling of CFT to two dimensional quantum gravity. 

The partition function of many quantum systems at in- 
verse temperature /3 may be rewritten as a Feynman path 
integral in imaginary time < r < /3. In the case of a cZ 
dimensional quantum lattice system we may regard the 
quantum partition function Z — ii e^^ as the canonical 
partition function of a classical system in d -|- 1 dimen- 
sions in a slab geometry with thickness (3. The boundary 
conditions in the imaginary time direction are periodic 
for bosonic degrees of freedom. The case when the clas- 
sical system is infinite in all d -f 1 directions corresponds 
to zero temperature in the quantum system. For sake of 
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FIG. 1: The critical ratio A"^ = C{XL, 8) /C{L, 8) as a func- 
tion of the fugacity z for L = 96, A = | and two rephcas (see 
eq.s (O and (H))). The top and the bottom arrows indicate 



respectively the expected value A 



for the pure system 



{z = 0) and A ^ for the coupled system at 2: = 1. 

and Aa are related through the KPZ formula. 



simplicity we assume that the classical system is defined 
on a d + 1 dimensional hyper-cubic lattice A = {x, r} 
{xi,T S Z). Its partition function can be computed by 
doing the Euclidean functional integral Z — J 'D[(j)]e^"^^'^^ 
over fields (j)x = 4>{x, r) periodic under t ^ t + j3. The 
Euclidean action S[(t>\ is assumed to be decomposable as 
the sum S = X]{a;j/) ^['t'x^ 't'y] of the contributions of the 
links ixy) of A . 

Working in the framework of replica method, we con- 
sider a stack of n copies of the original system, range 
them in a cyclical order and couple them together in the 
following way. We pick all the lattice nodes belonging to 
the inaccessible subsystem B and replace their links in 
the imaginary time direction with links connecting two 
consecutive copies 



c(fc) 



X e B 

x<^B 



(2) 



where <^'^^^ denotes the field in the fc— th copy. The sub- 
system B lies in a slice of the lattice at a given value of 
r. Since the imaginary time and space directions enter 
into the problem on very different footings, one might 
expect the corresponding classical system to exhibit in- 
trinsically anisotropic scaling. Here we restrict our anal- 
ysis to systems which are sufficiently isotropic, in such 
a way that their critical behavior is described by a rel- 
ativistic d + 1 dimensional field theory. In this context 
it is convenient to define a slightly generalized coupling 
among the n replicas, getting rid of the constraint of the 
B subsystem to lie in a constant slice and treating spa- 
tial and temporal links in the same way. Each stack of 
n links associated to the nodes x and y in the n replicas 



is set in two possible states. In the state 'A' each link 
of the stack connects points of the same replica while in 
the state 'B' it connects them cyclically, like in ([2|); the 
links in the state 'B' single out the subsystem B. An 
advantage of this more general setting is that it is easy 
to show that the coupled system of n replicas is endowed 
with an important local symmetry: flipping from 'A' to 
'B' or vice versa the state of all links intersecting an arbi- 
trary closed d dimensional manifold keeps invariant the 
partition function Q. A direct consequence of such a 
symmetry is that not only the entropy but all the ther- 
modynamic functions depend only on the boundary of B 
i- 

In order to promote accessible subsystems to dynami- 
cal variables one has simply to sum over all possible as- 
signments of the states 'A' and 'B' to the lattice links, so 
the partition function of our coupled system of n replicas 
can be written as 



Zn — 



„(fc) 



(3) 



where G is the subgraph of links which are set in the 
state 'B' and the summation is over all subgraphs. 

In the spirit of replica method, the evaluation of the en- 
tanglement entropy would require taking the limit n — ^ 1 , 
even if there are indications that the analytic continua- 
tion from positive integer n to real values could be rather 
difficult [S]. Fortunately, we do not need to do that, be- 
cause the new phenomenon we want to describe can be 
observed for any integer n > 1. 

To make the discussion concrete and explicit, we spe- 
cialize now to the case where the system in question is 
a quantum spin chain described at criticality by a CFT. 
Its two dimensional lattice description ([3]) is a discretized 
version of a n— sheeted covering of the plane, where the 
dual G of the subgraph G is formed by the set of cuts 
connecting these sheets. The local symmetry mentioned 
above turns out to express the invariance of the system 
under the addition (or the removal) of closed cuts or un- 
der continuous deformations of open cuts with fixed ends. 
It is worth noting that the only elements of the graph G 
having an intrinsic geometrical -and physical- meaning 
are the end points of the cuts, i. e. the branch points 
of the Riemann surface. They correspond to conical sin- 
gularities with deficit angle 27r(n — 1). On the contrary, 
the cuts joining different branch points are in no way 
distinguished lines on the surface: their introduction has 
a similar role as the choice of a reference frame on the 
surface. Thus, the dynamical effects of the back-reaction 
of the accessible subsystems of a critical quantum spin 
chain are intimately related to a 2D CFT in statistical 
equilibrium with a gas of conical singularities. 

Note that the n-sheeted covering of the plane with 
N branch points is a Riemann surface of genus g = 
{n — 1)(A^ — 2)/2, therefore summing over all accessi- 
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ble subsystems corresponds, in the replica approach, to 
a double sum over genera and moduli of these Riemann 
surfaces. This is the first indication that this issue is 
related to 2D quantum gravity. 

A stronger indication comes from an observation made 
by Knizhnik long time ago , namely that branch points 
correspond to primary fields $„(z,z) of scaling dimen- 
sions 

A«-A„^^(l--i) , (4) 

where c is the central charge of the system at the critical 
point. We see that for c not too large the conical sin- 
gularities are associated to relevant operators. Therefore 
the effect of considering the accessible subsystems as dy- 
namical quantities is equivalent, in the underlying field 
theory, to perturbing the CFT by the relevant operator 
$„(z, z), so the action is 

S = S* J ^n{z,z)(fz (5) 

where S* is the CFT action, and ^ is the chemical poten- 
tial which controls the appearance of conical singularities. 
Such a perturbation drives the system away form the 
critical point. According to the c— theorem [8[, a generic 
relevant operator generates renormalization group fiows 
into fixed points with a value of c that cannot exceed its 
initial value. However, there is more information to be 
gained. According to (H)) any CFT with c > is unstable 
whenever perturbed with This restricts the possible 
fixed points to which any CFT with c > may flow to 
those with c = 0. For generic values of these are trivial 
fixed points corresponding to massive theories. There is 
however one value of /x, at least, at which the system is 
critical. Thus the total central charge of the CFT cou- 
pled to the ensemble of the accessible subsystems is zero, 
which is precisely what it happens in CFT's coupled to 
quantum gravity. 

The strongest indication that this new universality 
class corresponds to 2D quantum gravity comes from nu- 
merical experiments, where one may accurately evalu- 
ate how change the scaling dimensions of local operators 
when the coupling to accessible subsystems is switched 
on. We use as a guidance the formula of Knizhnik, 
Polyakov and Zamolodchikov (KPZ)[i[i3l 

i".A4A(i-l), (6, 

which relates the scaling dimensions A° of a primary field 
of a CFT to the scaling dimensions A of the same oper- 
ator when the theory is coupled to quantum gravity. 

Let us consider a spin-^ quantum chain coupled to a 
transverse magnetic field h. The quantum Hamiltonian 
\s H = —J^iCr^af_^_i — h^-af , where and are 
the usual Pauli matrices. As is well known, this system 
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FIG. 2: The two bottom lines are one-parameter fits of 
(C(L,4) + C(L,8))/2 in the ID setting at z = 1 for n = 2 
and n = 3. The top line is a similar fit in the 2D setting at 
z — Zc for n = 2. In this doubly logarithmic plot the data lie 
on straight lines. The slope of the top line is exactly — 4Act 
while in the ID settings is shifted by the integer 1 — n. In 
order to represent all data in a single plot, the 2D data are 
reduced by a factor of 10. 

exhibits a quantum phase transition for J = h. This 
manifests itself as a power law decay of the correlators 
which lies in the universality class of the 2D critical Ising 
model, described by the CFT with c= ^. 

We simulated this system at the self-dual point with 
n = 2, . . . , 5 replicas of a square lattice enclosed in a 
square box with a side of L lattice spacings with toroidal 
boundary conditions. In a first set of Monte Carlo cal- 
culations the ensemble {A} was taken on a ID slice at a 
fixed value of r. On the same temporal slice we measured 
the (difference of) correlators 

= , (7) 

where s is a (proper) factor of the integer L. In our simu- 
lations we chose s = 8 and s = 4. In order to control the 
coupling with the ensemble {A}, we modified the Ising 
model by allowing a non-zero fugacity z — e^^ counting 
the number of branch points. 

Since {A} affected only a one dimensional boundary 
of our system, this remained critical even at z ^ 0. We 
simulated the Ising part of the model with a standard 
non-local cluster algorithm, while the update of the sub- 
set A, consisting of an arbitrary number N of disjoint 
intervals with N — 0, 1, ... , L/2, was performed with a 
heat bath method. Although the specific form of C{L, s) 
is not known, we have, in the thermodynamic limit, 

C{XL,s) = X-''C{L,s) (8) 

where A is any positive rescaling factor. As z varies from 
to 1, a; drops from the expected value of the pure Ising 
model X = 4A° to a new scaling dimension (see figure 1). 
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FIG. 3: The link variable as a function of the size of the 
lattice. The solid curve is a fit to eq.®. 



According to the previous discussion, the expected value 
is the one suggested by KPZ formula. As sometimes it 
happens in critical systems, the observed value of x differs 
by an integer with respect the expected value. Precisely 
we found x — 4Acr + 1 for n = 2 and x = AAr, + 2 for 



n > 3 (see figures 1 and 2). A^, 



is the KPZ value 



associated to A° = for 

lb 



Let us emphasize that in the first set of numerical ex- 
periments the simulated system is one dimensional, in the 
sense that the correlator and the subsystems are taken 
on a fixed ID slice at a given value of r. In the second 
set of numerical experiments we simulated instead a truly 
two dimensional system, with no limitations on the loca- 
tion of cuts representing the accessible subsystems. In 
this new setting the gas of conical singularities is spread 
in the bulk and drives the system away from the criti- 
cal point of the pure system. There is however a critical 
value Zc of the fugacity at which the whole interacting 
system undergoes a second order phase transition which 
is presumably in the same universality class of the ID 
quantum system described above. The system composed 
of two replicas at the self-dual point of the pure Ising 
system exhibits a critical behavior for Zc = 0.01127(1). 
The scaling dimension of the spin operator turns out to 
be, within the numerical accuracy, the one of the 2D 
quantum gravity (see figure 2). Likewise, the scaling di- 
mensions of the energy operator turn out to be the KPZ 
value Aj = |, corresponding to A° = i. In order to 
extract this critical exponent, we measured the vacuum 
expectation value of the link operator (cr| (jy), where x 
and y represent two nearest neighbors of the lattice. At 
the critical point this quantity is expected to have the 
following functional form 



where we wrote explicitly only the terms necessary 



(9) 



accurately fit the data, as shown in figure 3. 

From a geometrical point of view, the two-sheeted sys- 
tem we described is defined on an hyperelliptic Riemann 
surface and the sum X]{g} over all possible subsets of 
links made in eq.Q corresponds, in the continuum limit, 
to a double sum over the space of moduli and the gen- 
era of these surfaces. This has some similarity with the 
double scaling limit of matrix models llj. A promising 
aspect of the present approach is that it may be extended 
in a straightforward way to higher dimensions. 

Let us conclude with a general remark. In the pas t 
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years, 't Hooft [1^ and several other authors @, 
15| have suggested that the quantum entanglement and 
its ensuing entropy might be related to the Bekenstein- 
Hawking entropy of black holes. More recently, in the 
light of the AdS/CFT correspondence, a comprehensive 
gravitational interpretation of the entanglement entropy 
has been proposed The present study captures a 

different facet of the same fascinating relationship be- 
tween gravity and quantum entanglement, in which the 
back-reaction of the accessible subsystems of a ID quan- 
tum system has the same effect as the coupling to 2D 
quantum gravity. 
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